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Abstract 

An aerodynamic database has been generated for 
the Mars Smart Lander Shelf-All configuration using 
computational fluid dynamics (CFD) simulations. Three 
different CFD codes, USM3D and FELISA, based on 
unstructured grid technology and LAURA, an estab- 
lished and validated structured CFD code, were used. 
As part of this database development, the results for the 
Mars continuum were validated with experimental data 
and comparisons made where applicable. The valida- 
tion of USM3D and LAURA with the Unitary experi- 
mental data, the use of intermediate LAURA check 
analyses, as well as the validation of FELISA with the 
Mach 6 CF 4 experimental data provided a higher confi- 
dence in the ability for CFD to provide aerodynamic 
data in order to determine the static trim characteristics 
for longitudinal stability. The analyses of the non- 
continuum regime showed the existence of multiple 
trim angles of attack that can be unstable or stable trim 
points. This information is needed to design guidance 
controller throughout the trajectory. 
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Nomenclature 


A 

Reference area (m 2 ) 

BC 

Ballistic coefficient (m /CD*A) 

CA 

Axial force coefficient 

CD 

Drag coefficient 

CG 

Center of gravity 

CN 

Normal force coefficient 

Cm 

Pitching moment coefficient 

Kn 

Knudsen number 

LID 

Lift to drag ratio 

m 

Vehicle mass (kg) 

M 

Mach number 

a 

Angle of attack 

V 

Ratio of specific heats 

Cm a 

Static stability derivative 

Subscripts: 


trim 

trim condition 



Introduction 

The exploration of Mars has been driven by the 
need to understand its planetary evolution and thereby 
recognize its past, present, or future potential to support 
life. The discovery of Martian life would require the 
ability to understand or characterize Mars from its ge- 
ology, climate, biology, and other natural processes. 
This requirement ultimately defines a set of overall 
science objectives/payloads along with a probable set of 
planetary surface locations. In addition, the accuracy of 
targeting a specific surface location and achieving a 
successful precision landing are also essential require- 
ments for fulfilling these science objectives. 

The entry, descent, and landing (EDL) phases for a 
Mars mission can be categorized based on 1 st , 2 nd , and 
3 rd generation type systems. 1 These categories describe 
the development of landing ellipses based on large scale 
(100s of km), precision (< 10 km), and pinpoint dis- 
tance accuracies from the desired landing site. The I s1 
generation EDL system utilizes a ballistic entry. The 2 nd 
generation EDL system (Smart) employs significant 
improvements in guidance and navigation methods 2 ' 3 
enabling precision entry and utilizes a “Smart” or local 
surface hazard avoidance and hazard tolerance system 
in order to provide a safer landing environment. The 2 ,ld 
generation EDL system also requires the development 
of a higher fidelity trajectory simulation. The objectives 
are for the vehicle to land within 10 km of the intended 
touchdown site and landing safely 99.7% of the time. 

The trajectory simulation for the Mars Smart 
Lander (MSL) configuration series was performed us- 
ing the 3-degree of freedom (DOF) version of the 
POST 4 (Program to Optimize Simulated Trajectories). 
The POST program integrates the equations of motion 
and includes a set of specific constraints necessary to 
define the overall planetary, mission, and vehicle pa- 
rameters. These models include specific science-related 
information of the planet Mars such as planet defini- 
tion, gravity, and atmosphere while other models pro- 
vide specific information related to the vehicle configu- 
ration and the vehicle characteristics of motion. Exam- 
ples of these models include control system emulation, 
guidance algorithms, navigation system, mass proper- 
ties, and aerodynamics. The aerodynamics model is a 
database that provides aerodynamic force and moment 
data of a specific vehicle configuration for both the 
entry and descent segments of the trajectory. The entry 
phase is defined as the region between the atmospheric 
interface and the supersonic parachute deployment 
point while the descent phase is defined from that point 
to landing. 


The development of an entry aerodynamic database 
requires specific knowledge of the geometry and ori- 
entation of the vehicle in addition to the freestream 
quantities along a given trajectory. In accordance with 
the Mars Smart Lander payload and surface location 
requirements, a set of preliminary configurations were 
developed that would have favorable trim characteris- 
tics along proposed nominal and dispersed trajectories. 
These vehicle configurations (Figure 1) are the axi- 
symmetric Baseline, Canted-All, Shelf-All, and a de- 
rivative of the Shelf-All tested in the Ames Ballistic 
range facility. All configurations have the same 70° 
forebody cone shape with a bi-conic backshell. The 
Shelf-All model was designed with a blended control 
tab that is tangent to the forebody surface while the 
Shelf-2 model has a non-blended control tab while also 
being tangent to the forebody surface. The Canted-All 
model has a similar Shelf-All control tab but is canted 
10° towards the vehicle nose. 



Figure 1. Mars Smart Lander Configurations 

The objectives of this paper are to provide infor- 
mation on the development of an entry aerodynamic 
database, to provide specific validation of the CFD 
analysis programs with experimental data, and show the 
resulting flight trim characteristics for a Mars Smart 
Lander configuration. The subsequent sections of this 
paper will describe the analysis codes used, CFD vali- 
dation, and the creation and results of a flight database. 

Code Descriptions 

USM3D 

USM3Dns 5 is a three dimensional, tetrahedral, 
cell-centered, finite volume Euler and Navier-Stokes 
flow solver for unstructured meshes. Inviscid flux 
quantities are computed across each cell face using 
Roe’s flux-difference splitting (FDS) 6 Spatial discreti- 
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zation is accomplished using an analytical reconstruc- 
tion process for computing solution gradients within 
tetrahedral cells. The solution is advanced in time to a 
steady state condition by an implicit backward Euler 
time-stepping scheme. Flow turbulence is modeled by 
the Spalart-Allmaras one-equation model, which is op- 
tionally coupled with a wall-function formulation to 
reduce solution stiffness and the number of cells in the 
sub layer of the boundary layer. USM3Dns runs on 
massively parallel computers and clusters of personal 
computers (PC’s). Although a single processor version 
is available for a variety of computing platforms, the 
parallel version 1 2 3 * * * 7 is the code of choice because it en- 
ables rapid turnaround for large problems. 

Although USM3Dns is widely in use by aerospace 
practitioners for several years, the latest Navier-Stokes 
version was never validated for the Mach number range 
and shape (especially the blunt fore body) of interest in 
the present study. In addition, the available version was 
so far used with perfect gas assumptions, while the 
aerodynamic database required solutions at Mars at- 
mospheric conditions involving a different gas constant. 
For these reasons, the flow solver had to be modified 
before application. This section briefly describes the 
changes made to the code. The following changes were 
incorporated to the existing parallel version of 
USM3Dns: 

1. A new flow initialization procedure was intro- 
duced whereby for the cases were the free- 
stream Mach number exceeded 1.2; the solution 
was initialized to a lower, usually subsonic, 
value. This strategy helped the flow solution 
convergence to steady state, especially in the 
aft, low-Mach number region of the configura- 
tions. 

2. A special flux-splitting scheme based on the 
AUSM 8 * * * (Advection Upstream Splitting 
Method) scheme was incorporated to overcome 
the well-known “Carbuncle phenomenon” usu- 
ally associated with the application of the Roe 
scheme to a case with strong shock waves. 

3. In order to simulate the Mars atmospheric con- 

ditions, the perfect gas assumptions used thus 

far was modified to allow specification of an ef- 

fective y for each freestream Mach number, 

along with coefficients for Sutherland’s for vis- 

cosity. 

LAURA 

The Langley Aerothermodynamic Upwind Relaxa- 

tion Algorithm (LAURA) is a three dimensional Na- 


vier-Stokes solver developed by Gnoffo 941 for struc- 
tured grids. LAURA is a finite volume formulation of 
the integral form of the Navier-Stokes equations. A 
second order-accurate, symmetric total-variation- 
diminishing (TVD) scheme 12 is used in conjunction 
with upwind differencing of the discretized equations. 
At each cell face. Roe’s averaging defines the flowfield 
variables based on values from the adjacent cells. The 
unsteady governing equations are driven to a steady 
state solution through an implicit time relaxation proce- 
dure. During the relaxation process, the grid is periodi- 
cally adapted in the body-normal direction so that the 
grid can be tailored to the emerging solution. This el- 
liptic flow solver has the capability of solving chemical 
non-equilibrium flow through the specification of 
chemical species and reactions unique to a specific at- 
mosphere as well as the ability to model perfect gas and 
equilibrium air. 

FELISA 

The FELISA 1344 (Finite Element Langley Imperial 
College Swansea Ames) software system consists of 
unstructured surface and volume grid generation with 
an inviscid flow solver and unstructured post process- 
ing utilities. The flow algorithm was developed through 
the application of the Galerkin finite element method in 
space to obtain a coupled set of ordinary differential 
equations in time. The steady state solution of this 
equation set is achieved by advancing the system using 
an explicit Runge-Kutta type marching scheme. The 
hypersonic flow solver has options for perfect gas, 
equilibrium air, CF 4 , C0 2 , and Mars equilibrium gases. 
In addition, the hypersonic flow solver has the capabil- 
ity of solving chemical non-equilibrium flow and real 
gas (chemical and thermal non -equilibrium) flow. 

Aerodynamic Database 

The aerodynamic database is provided as a 
FORTRAN routine containing specific vehicle force 
and moment data as a function of vehicle orientation 
relative to the flow and flowfield parameters to the 
POST trajectory simulation. Within this routine, the 
vehicle trajectory is partitioned into three atmospheric 
regimes that are defined based on the similarity pa- 
rameter of Knudsen number. The Knudsen number 
(Kn) is defined as the mean free path (distance between 
molecular collisions) divided by the vehicle reference 
length (typically the aeroshell diameter). These atmos- 
pheric regimes are depicted in Figure 2. 
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Smart Lander Reference Trajectory 
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Figure 2. Mars Smart Lander Reference Trajectory 


The aerodynamic database is comprised of contin- 
uum, transitional, and free molecular 15 16 regimes that 
are accessed based on the calculated local Knudsen 
number at any entry trajectory point. Each section con- 
sists of axial force, normal force, and pitching moment 
data that was created from specific CFD analysis points 
as well as experimental data 1 that have been corrected 
for testing environment differences. Since the trajectory 
simulation is normally run under nominal and off- 
nominal conditions during prescribed Monte Carlo 
analyses, the aerodynamic database has to be suffi- 
ciently populated at several freestream Mach numbers 
and angles of attack. To determine values within or 
outside the Mach number and angle of attack ranges, a 
bi-parabolic/linear interpolation scheme is applied 
within the range of the database. 

Within the continuum regime, three different CFD 
analysis programs were used to obtain direct flight 
aerodynamic data or to provide necessary corrections to 
the Unitary experimental data for flight condition appli- 
cability while the non-continuum flight data was deter- 
mined using the program as mentioned in References 
15 and 16. The USM3Dns code was used to perform 
the Unitary data corrections from the predicted low 
supersonic parachute deploy boundary to the upper su- 
personic limit of the Unitary facility. The FEL1SA re- 
sults provided the aerodynamic flight results from a 
mid-hypersonic to upper hypersonic trajectory point 
that corresponded to the approximate atmospheric con- 
tinuum limit. Lastly, the LAURA results provided in- 
termediate analysis comparisons that overlap the above 
trajectory ranges. The proposed continuum flight cases 
are shown in Table 1. 


The selections and application of these analysis 
programs were based on several considerations. The 
first consideration was to develop a methodology that 
could allow for rapid turnaround of assessing trends for 
proposed vehicle configurations by taking advantage of 
grid generation methods and flowfield modeling. An- 
other consideration was the flight regime capability of 
each program. The final consideration was to ensure 
confidence and continuity within the solution matrix by 
using the LAURA code to validate intermediate analy- 
sis points with the USM3D and FELISA analysis pro- 
grams since LAURA has been previously validated 
with previous Mars missions 18 . 


Table 1. Proposed Flight Analysis Cases for MSL 
Reference Trajectory. 

Analysis Program 

M 

V (m/s) 




1.6 

361.37 

USM3D 



2.3 

510.38 



3.5 

755.52 




4.5 

957.06 


LAURA 


6 

1262.97 



8 

1677.90 




10 

2092.29 




13 

2706.78 



FELISA 

16 

3298.42 



20 

4039.48 




24.5 

4763.74 




28 

5295.34 


Experimental Data 

Two experimental facilities were used to acquire 
aerodynamic data for CFD validation of the aforemen- 
tioned analysis programs. The first facility, the Langley 
Unitary Plan Wind Tunnel, is a closed-circuit, continu- 
ous flow, variable density supersonic tunnel that was 
used to obtain 6-DOF static aerodynamic data for 4 
proposed Mars Lander configurations over a Mach 
number range of 2.3 to 4.5. The complete details of the 
models, facility calibration, instrumentation, and data 
acquisition are presented in Reference 17. The second 
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facility, the 20-inch Mach 6 CF 4 Tunnel, was used to 
obtain aerodynamic data for the Mars Surveyor con- 
figuration. The details of this work are presented in 
Reference 19. 


Results and Discussion 

USM3Dns Validation 

The USM3D code was required to validate its re- 
sults with Unitary tunnel experimental data in order to 
provide confidence in correcting the Unitary Wind tun- 
nel data for the geometric differences between wind 
tunnel and flight configuration and gas composition of 
the Mars atmosphere. Along with the previously de- 
scribed changes to USM3D, an unstructured mesh was 
constructed for the Baseline and Shelf-All tunnel con- 
figurations for a series of Mach numbers between 2.3 
and 4.5, each for a range of angles-of-attack values. The 
grid for the full Navier-Stokes solution varied in size 
typically from 1.4 million tetrahedral cells for the wind- 
tunnel cases to 2.6 million cells for the flight configu- 
rations. All the calculations were run for a Reynolds 
number of 1.0 million based on the maximum diameter 
of the configuration (6 inches for the wind tunnel mod- 
els and 4.05 meters for the flight cases). The CFD cases 
were run using up to 60 processors on an SGI Origin 
2000 computer. For all the computations reported here, 
a full viscous formulation was used. Some examples of 
the unstructured surface topology are depicted in Fig- 
ures 3 and 4. 



Figure 4. Unstructured Surface Grid for the Unitary 
Shelf- All Model Configuration. 

Figure 5 for the Baseline configuration shows the 
USM3D/Unitary comparisons for freestream Mach 
numbers of 2.3 and 3.5. The axial and normal forces 
along with the pitching moment (about the nose) are 
shown. The CFD results for the tunnel conditions (y=1.4) 
compare well against the Unitary wind tunnel data, at- 
testing to the accuracy of the modified solver. The de- 
tailed comparisons show a better overall agreement at 
the Mach 3.5 condition then at the lower Mach 2.3 con- 
dition while the Mach 2.3 condition also depicting 
smaller differences in axial and normal force at 0 and 20 
degrees than the intermediate angle of attack range. 



Figure 3. Unstructured Surface Grid for the Unitary 
Baseline Model Configuration. 
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Figure 5. USM3D/Unitary Data Comparisons for the 
Baseline Configuration. 


However, the USM3D analysis does predict a change in 
normal force slope near zero degrees angle of attack 
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which is a known trend for axisymmetric blunt body 
configurations at low supersonic speeds. 

Figure 6 shows the comparisons for the Shelf-All 
configuration to the Unitary experimental data at 
freestream Mach numbers of 2.3, 3.5 and 4.5. As with 
the previous comparisons shown for the Baseline con- 
figuration, the overall agreement for Shelf-All configu- 
ration is very good. These comparisons show better 
agreement in axial and normal force at lower Mach 
numbers than shown for the Baseline model while the 
pitching moment and L/D distribution shows excellent 
agreement throughout the Mach number and angle of 
attack range. 
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Figure 6. USM3D/Unitary Data Comparisons for the 
Shelf-All Configuration. 

Figure 7 shows the comparisons for the Shelf-All 
configuration to the Unitary experimental data under a 
sideslip condition of 4 degrees. As with the previous 
comparisons shown for the Baseline configuration, the 
overall agreement for Shelf-All configuration is very 
good. This particular case (6-DOF) was chosen in order 
to increase the level of validation and thereby illustrate 
any deficiencies due to a geometry change and under 
sideslip conditions. These comparisons show a better 
agreement in axial force and a similar difference in 
normal force as compared to the non-sideslip 
freestream condition. Likewise, the comparisons to the 
pitching moment and L/D are excellent. Further as- 
sessment of CFD validation would require some solu- 
tions at lower angles of attack, but the overall assess- 
ment is extremely good. It should also be noted that this 
sideslip case was not used in the creation of the flight 
database that is presented in this study. 
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Figure 7. USM3D/Unitary Data Sideslip Compari- 
sons for the Shelf- All Configuration (P=4°). 


Figure 8 shows the comparisons for the Canted-All 
configuration to the LInitary experimental data at 
freestream Mach numbers of 2.3 and 3.5. As with the 
previous comparisons shown for the Baseline configu- 
ration, the overall agreement for Canted-All configura- 
tion is very good. These comparisons show better 
agreement throughout the Mach number and angle of 
attack range for all force, moment and L/D quantities 



Figure 8. USM3D/U nitary Data Comparisons for the 
Canted-All Configuration. 


than those shown for the Baseline configuration. FIow- 
ever, additional solutions at lower angles of attack 
would have to be made in order to obtain a more com- 
plete assessment of CFD validation. 
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LAURA Validation 


Although the LAURA analysis program has been 
previously validated with Mars Pathfinder hypersonic 
flight data 18 , an attempt was made to provide additional 
validation with the Unitary experimental data for a sin- 
gle Mars Smart Lander configuration at supersonic 
speeds and to further validate USM3D. Figure 9 shows 
the LAURA structured grid topology for the Unitary 
Baseline configuration. This structured grid is com- 
prised of 16 computational blocks with an approximate 
total of 450,000 cell centers. In addition, this structured 
grid was adapted in the direction perpendicular to the 
body surface as a function of the flowfield temperature 
and grid distribution values. 




Figure 10. Mach Number Contours for the Unitary 
Baseline Configuration - (M=4.5, a=10°) 


Figure 11 depicts the LAURA/Unitary compari- 
sons for the Baseline Unitary tunnel configuration for 
freestream Mach numbers of 2.3 and 4.5. These com- 
parisons overall show good agreement to the Unitary 
data for all quantities. The axial force comparisons 
show better agreement at the higher angles of attack for 
both Mach numbers while the normal force comparison 
show-s better agreement over the entire angle of attack 
range at Mach 4.5 while the 10 degree case for Mach 
2.3 show's a larger difference as compared to the 
USM3D/Unitary comparisons in Figure 5. 


Figure 9. Structured Grid for the Unitary Baseline 
Model Configuration 


Figure 10 depicts the LAURA Mach number con- 
tours for the Baseline Unitary tunnel configuration for 
freestream Mach number of 4.5 and an angle of attack 
of 10 degrees with no sideslip angle. This solution im- 
posed perfect gas air (y=1.4) at M=4.5 tunnel conditions 
with a laminar, no-slip adiabatic wall boundary condi- 
tion. Figure 10 shows along the vehicle symmetry and 
outflow planes the general flowfield characteristics of 
an axisymmetric blunt body with dominant subsonic 
regions within the shock layer and vehicle base. 



Figure 11. LAURA/Unitary Data Comparisons for 
the Baseline Configuration. 
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FELISA Validation 


The FELISA validation work was centered about 
an earlier Mars 2001 Surveyor configuration that was 
tested in the 20-inch Mach 6 CF 4 Tunnel. The Mars 
Surveyor configuration (Figure 12) was similar to Mars 
Smart Lander configuration, but differed in the control 
flap geometry and location. The Mars Surveyor flap 
surface area was at least 50 percent smaller and the flap 
orientation was 90 degrees perpendicular from the axis 
of rotation and located at aft of the maximum diameter. 
The number assigned to each control flap indicated the 
area ratio of the flap where the larger number has a 
greater area ratio. 



Figure 12. Mars 2001 Surveyor Wind Tunnel Model 

Figures 13 and 14 depict the pitching moment co- 
efficient and L/D ratio for the Mars Surveyor configu- 
ration for Flap 1. The FELISA CF 4 results showed very 
good agreement to the experimental CF 4 data compari- 
sons at 0, 5, 11, and 16 degrees angle of attack. In con- 
junction with a planned test in the 20-inch Mach 6 fa- 
cility, a FELISA analysis was performed using a C0 2 
test gas as a freestream condition. The FELISA results 
predicted that the force and moment contributions to the 
L/D ratio at 0 degrees angle of attack would be insensi- 
tive to the C0 2 test gas. However, the need for addi- 
tional results at different angles of attack would be nec- 
essary in order to determine a more generalized trend in 
C02 insensitivity. 



Figure 13. FELISA/CF 4 Pitching Moment Compari- 
sons for the Mars Surveyor Configuration. 



Figure 14. FELISA/CIj L/D Comparisons for the 
Mars Surveyor Configuration. 


Figures 15 and 16 depict the pitching moment co- 
efficient and L/D ratio for the Mars Surveyor configu- 
ration with Flap 3. The FELISA CF 4 results shows very 
good agreement to the experimental CF 4 data compari- 
sons at 0, 5, 11, and 16 degrees angle of attack. The 
larger Flap 3 configuration does show a slightly larger 
deviation in pitching moment at 0 and 5 degrees as 
compared to the Flapl, but overall agreement is still 
very good. These differences could be attributed to vis- 
cous contributions (not being modeled in FELISA) in- 
duced by the larger flap or an indication of a required 
increase in control flap grid resolution. 
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determination of ballistic coefficient, a mass of 2200 
kg. was used. 


C,. Vs oc FOR '01 MARS LANDER FOREBODY IN CF, 
Flap 3 



At the axial CG of 0.21, the trim angle of attack 
depicts oscillatory behavior of 1 degree within the 
Mach 2.3 to 4.5 LInitary experimental data range before 
reaching a maximum value of 15.4 degrees at a velocity 
of approximately 1300 m/s (Mach 6). Beyond Mach 6, 
the trim angle of attack at the forward CG decreases 
nearly linear to a continuum limit value of approxi- 
mately 12.8 degrees. Likewise, the L/D values depict 
the same low supersonic behavior as well as the nearly 
linear behavior beyond the maximum L/D of 0.266 at 
Mach 4.5. 


Figure 15. FELISA/CF 4 Pitching Moment Compari- 
sons for the Mars Surveyor Configuration. 



Figure 16. FELISA/CFJ L/D Comparisons for the 
Mars Surveyor Configuration. 


Flight Results - (Shelf- All Continuum Database) 

For the continuum regime, an aerodynamic data- 
base was created for the Mars Smart Lander “Shelf-All” 
configuration using the USM3Dns correction analyses 
for the Unitary experimental data (M=2. 3-4.5), the 
LAURA intermediate analysis check points at Mach 4.5 
and 10, and the FELISA analyses at Mach 6, 8, 10, 24, 
and 28. Figures 17 and 18 show the continuum trim 
characteristics of the Shelf-All configuration based on 
the axial CG locations of (x/Diameter) of 0.21 and 0.32. 
The CG of 0.21 corresponds to an early estimate of the 
MSL CG location while the 0.32 value corresponds to 
an L/D value of approximately 0.24 at a velocity of 
4700 m/s (Mach 24). In addition, the radial CG offset 
for this configuration was designed to be zero. For the 


At the axial CG of 0.32, the trim angle of attack 
depicts a larger variation in trim angle of attack of 
about 2.5 degrees within the LInitary data range and 
also reaching a maximum of 16.9 degrees at a velocity 
of approximately 550 m/s (Mach 2.7) while reaching 
another lower peak of 16.1 degrees at Mach 6. Beyond 
Mach 6, the trim angle of attack at the 0.32 CG de- 
creases nearly linear to a continuum limit value of ap- 
proximately 13.6 degrees. As with the forward CG re- 
sults, the L/D shows a corresponding maximum value 
at the peak trim reaching an overall maximum of 0.298. 
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Figure 17. L/D Continuum Trim Characteristics for 
Shelf- All Configuration. 


At trim conditions, the ballistic coefficient for the 
forward CG shows no low supersonic oscillatory be- 
havior while approaching a maximum of 1 14 at Mach 6 
while the axial CG of 0.32 showed the same distribu- 
tion while displaying a maximum of 116 at Mach 6. As 
expected, these results show a low sensitivity to CG 
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location due to only a small (1%) difference in ballistic 
coefficient over the entire continuum regime. 



Figure 18. Ballistic Coefficient Continuum Trim 
Characteristics for Shelf- All Configuration. 


Figure 19. Non-Continuum Pitching Moment for 
Shelf- All Configuration. 



Flight Results-(Shelf-All Non-Continuum Database) 

For the rarefied regime, an existing set of free mo- 
lecular data 2 over an angle of attack range of -180 to 
180 degrees was used in conjunction with bridging 
function 20 to provide aerodynamic values within the 
transition atmospheric regime. Figures 19 and 20 show 
the trim characteristics for the Shelf-All configuration 
from the edge of continuum to a value of 20. A Knud- 
sen number (Kn) of 20 was sufficient to cover the range 
up to the atmospheric interface at approximately 
125000 meters above ground level. 

Figure 19 shows the oscillatory behavior of the 
pitching moment within the transition regime. At the 
continuum limit and up to a Knudsen number of ap- 
proximately 5, there exists only a single trim point that 
is also stable (negative Cm a ). As the vehicle descends 
through the transition regime, there are multiple condi- 
tions that the vehicle could potentially trim, but could 
potentially be an unstable or stable trim point. Figure 20 
displays Knudsen number and Cm alpha as a function 
of trim angle of attack (a trim ). For Knudsen numbers 
above 2, there exist three trim angles of which 2 are 
unstable and 1 is stable. This information is needed 
while designing a guidance controller during the entry 
phase of the trajectory. 
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Concluding Remarks 

The CFD validation of USM3Dns and LAURA 
with the Unitary experimental data, the use of interme- 
diate LAURA check analyses, as well as the validation 
of FELISA with the Mach 6 CF 4 experimental data pro- 
vided very good agreement and a higher confidence in 
each CFD code’s aerodynamic analysis results includ- 
ing the longitudinal static trim characteristics for the 
Mars continuum atmosphere. 
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